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, Abstract. This is a review of formalisms and models (nonrelativistic and relativistic) 

which modify Schrodinger's equation so that it describes wavefunction collapse as a 
On ' dynamical physical process. 

(N 

, 1 Representing The World 

, One of the major goals of theoretical physics is to find a useful mathematical 

' object to represent the physical state of our world. In classical physics the object 

is a point in phase space. So far it has been possible to represent the world at an 
' "instant," i.e., on a spacelike hypersurface^. One thing "useful" means is that a 

dynamical equation can be found for the object, enabling predictions about the 
state of the world on a future spacelike hypersurface. 

Is there such an object in quantum physics? In pursuing the research dis- 
cussed here I have made some bets as to the nature of an eventually satisfactory 
physical theory. One of them is that there is such an object, a statevector in 
' a suitable Hilbert space plus something more. I shall argue that more must be 

^ , added because standard quantum theory (SQT) is a theory of choices without 

a chooser: more is a chooser. Along the way I shall make a few opinionated re- 
marks about other points of view. However, mostly I shall try to give a coherent 
account of rationale (section 1), formalisms (section 2), nonrelativisitic (includ- 
ing gravitational, section 3) and relativisitic (section 4) models for a particular 
' way of altering Schrodinger's equation so the statevector behaves in a manner 

consistent with this bet. 

There is a point of view that the statevector does not represent the world but 
is only a mathematical tool enabling people to calculate the statistical outcomes 
of experiments[l]. There is a point of view that only an ensemble of experiments 
is represented by the statevector [2, 3, 4]. I disagree with both because I find 
it hard to believe such a successful theory as quantum theory either has no 
relevance to the unobserved world or to the individual world. There is a point 

^ One can imagine other possibilities, e.g., a representation might only be possible 
on the spacetime volume between two hypersurfaces; there might be another time 
parameter whose specification is also required; etc. 
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of view that it is the density matrix which is the proper object to represent the 
world. There is a point of view that only ensembles of worlds can be represented. 
I shall comment on these later on in this section. 

I have of necessity been loose about what I mean by the part of "the world" 
whose state the statevector represents. It should contain at least what we observe 
exists and what we thereby infer exists. I will take as a minimal characterization 
that this means all particles and structures whose existence physicists at present 
generally agree upon. These do not include the "more" mentioned above, which 
is posited to be an additional structure in the world. However, it may be possible 
to include the "more" in the statevector as well (sec section 2.3h). 

Essential is the view that there is a specifiable "real world out there" [5] whose 
behavior is completely consistent with our observations. Implied is that we are 
fair observers (what we see is there) and sound reasoners about the world so 
our conclusions can be trusted. Not implied is that we are in any way essential 
to the description of the world. Rather it is the theory of the world that is 
essential to the description of us and our observations (what has been called 
empirical reality[6]). I bet against the ultimate value of any theory whose logical 
consistency requires human beings [7]. 

All statevectors having the same "direction" in Hilbert space characterize the 
same state of the world. This generalizes the "ray" (overall phase) independence 
of the SQT description to norm independence. In SQT the norm (squared) must 
bo 1 because it represents the sum of probabilities, but that is not required when 
the statevector represents actuality. If two statevectors differ in direction by the 
slightest bit they represent two different states of the world. However, any one 
of a "small" ball of directions will suffice to describe the same state of empirical 
reality [8, 9], i.e., an object may be empirically here when a small amount of its 
wavefunction is elsewhere. 

It is very possible that the statevector should represent more particles and 
structures than physicists now agree exist. However, I don't consider that the 
statevector should represent more variant replicas of our world, as in the many- 
worlds interpretation of the statevector[10, 11] because it is likely that the foun- 
dations of a coherent many- worlds picture requires human beings[12, 13, 14]. 
But, even with such an interpretation, one of the statevectors in the many- 
worlds superposition describes our world, and it is this statevector I want to 
talk about. 

The assumption that the statevector gives a complete specification of the fa- 
miliar particles and structures in the world excludes from discussion the deBroglie- 
Bohm model[15, 16, 17, 18, 19]. This model is akin to the models I discuss in 
that both contain an extra structure, a chooser, but the deBroglie-Bohm chooser 
is the positions of particles. I am betting that the statevector alone satisfactorily 
describes the particles in the world. Both the many-worlds and deBroglie-Bohm 
approaches to the central problem addressed here are interesting and have suc- 
cesses and difficulties, but I won't discuss them. 

This central problem is arguably the most important unsolved problem in 
the foundations of quantum mechanics. It is brought about by certain features 
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of the world as evidenced by our observations and inferences. These are that 
events occur in the world, that in certain circumstances they occur in a fun- 
damentally unpredictable (random) way, and that macroscopic objects (choose 
your own definition) arc almost always highly spatially localized — even objects 
whose location is an unpredictable event. 

This latter behavior means that certain statevectors cannot represent possible 
states of the world, namely statevectors whose evolution describes a macroscopic 
object continually in a superposition of here and there with comparably sized 
coefficients. The problem arises because, using the otherwise estimable SQT, 
one may start with a statevector representing a state of the world, in a situation 
such as a measurement in which uncontrollable events occur, and it evolves via 
Schrodinger's equation into a statevector which does not represent a state of the 
world, but rather the sum of such states. 

The founding fathers of quantum theory had a cure for this. It is to modify 
the evolution of the statevector so that it instantly "collapses" (or "jumps" or 
"reduces") in such circumstances to a viable statevector representing the state 
of the world. The problem with this modification is that it is terribly ill-defined. 
It is supposed to be invoked whenever a measurement has been completed, but 
no one has been able to define what a measurement is or when completion is. 
Ad hoc means "for this case only," and the prescription given is very ad hoc, 
with every situation requiring its own assessment. Perhaps due to the persuasive 
powers of John Bell [7], more physicists nowadays think this collapse postulate 
is unsatisfactory. 

There have been two approaches which can be said to follow in the footsteps 
of the founding fathers. 

One approach is to try to make the instantaneous collapse postulate well 
defined without the need for ad hoc information. This is my view of the "con- 
sistent histories" program[20, 21, 22]. It has been described as a promissory 
note[23] which has serious problems in being fulfilled[24]. The problem is that 
Nature selects a unique set of consistent collapse possibilities but the theory does 
not. 

There is a test I think should be appled to all theories with fundamental 
pretensions [25]. If, confronted with an initial statevector of a seriously compli- 
cated and realistic part of the world (e.g., the local galactic group), are there 
well-defined procedures for constructing the mathematical quantities which cor- 
respond to the possible real events which take place in the world? The point 
of this test is to see if a theory can do more than just handle some simplified 
models into which ad hoc features creep. So far the consistent histories approach 
does not pass this test. 

This is as good a place as any to mention that the somewhat related "en- 
vironmental decoherence" scheme[26, 27] does not pass this test either [28]: so 
far the choice of "environment" is ad hoc. This scheme is a prime example of 
trying to use the density matrix to represent the world, which has the the fol- 
lowing problems. If the density matrix is, in the Victorian characterization of 
d'Espagnat[29] "Pure and Proper", i.e., = \^,t >< tp,t\, then it is equivalent 
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to a statevector and must undergo collapse in order to represent the world. If 
it is Impure (mixed) and Proper then it may either represent an ensemble of 
worlds or our ignorance about our world, but not the world itself. But if it is 
e.g., a density matrix traced over the "environment," it begins as a Pure and 
Improper representation of our world and ends up interpreted as an Impure and 
Improper representation of an ensemble of worlds or of our ignorance about the 
world. The interpretation of a consistent theory should not undergo dynamical 
evolution. This is why I do not believe a density matrix is the appropriate object 
to represent the world. 

2 Describing Collapse 

A second approach following the founding fathers, and espoused here, is to agree 
with them that the evolution of the statevector should be modified but to do it 
so the collapse is not instantaneous but, rather, follows a well-defined dynamics 
of a modified Schrodinger equation. This is a bet that there is a real physi- 
cal process which causes events to occur which is not yet in physics and it is 
worthwhile to try to make a phenomenological model of it. Such a model is very 
strongly constrained because it must agree with empirical reality, which includes 
all of tested physics, the random choices made by nature and the localization of 
macroscopic objects. 

In order to describe the random choices made by nature the Schrodinger 
equation must have a chooser in it. The first such model[30] used some hidden 
variables. Subsequently, I tried to use the phases of the superposition[31] as 
chooser (which is appealing because nothing "more" must be added beyond 
what is in the statevector). Then, betting that the chooser is one of many things 
in nature which fluctuate randomly, I settled on modelling it by external random 
noise[32], a choice which has been adopted in subsequent work[33, 34, 35, 36, 37]. 
In this section I shall discuss variations on the theme of collapse formalism. 

2.1 Gambler's Ruin Game 

With external random noise as the chooser, it turns out that the mechanism for 
obtaining agreement with the predictions of quantum theory is very simple[38] 
which suggests one is on the right track. 

The mechanism is the gambler's ruin game[39]. Suppose one gambler has $36 
and another has $64, and they toss a coin to determine who gives a dollar to 
the other. Their dollar amounts fluctuate. Eventually one gambler wins all the 
money. It is readily shown that the probability of winning all the money is .36 
for the first (.64 for the second). 

Precisely analogous is the modelled evolution of the initial statevector 



IV',0 >= .6|ai > +.8|a2 > 



(1) 
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Under the influence of the random noise, the amplitudes multiplying |oi > and 
\a2 > fluctuate. Eventually the statevector ends up as |ai > with probability .36 
(ends up as \a2 > with probability .64). This final result is just what the SQT 
collapse postulate accomplishes. 

2.2 Simple Model 

Here is a simplest example I can give of such a collapse dynamics. 

Two equations are required for a collapse model. The first is the dynamical 
equation which replaces Schrodinger's equation. Remarkably, it is possible to 
write this as a linear equation for the statevector [40, 41]. This simplicity is 
very useful e.g., it makes possible the relativistic collapse models discussed in 
section 4. The dynamical equation depends upon a classical white noise function 
w(t). However, for this simple example, the solution actually depends only upon 
jQw{t')dt' = B{t) which is a classical Brownian motion function. Here is the 
solution \il),t >B when the initial state is (1): 

= .6|ai > e-jfet^W-^^*"*^!' + .8|a2 > e-^fet^W-^^*'*^]' . (3) 

As will soon be seen, the operator A determines the choices while B{t) is the 
chooser. A's eigenstates, |ai >, |a2 > (eigenvalues ai, 02) are the states to which 
collapse occurs. A's eigenvalue differences, together with A (a parameter of the 
theory) determine the collapse rate. 

The evolution equation (2) tells us what the initial statevector evolves into 
under a particular B{t). The evolution is not unitary, so statevectors evolving 
under different i3(t)'s have different norms. This plays a role in the second re- 
quired equation, the probability rule. It gives the probability density for B{t) to 
be the actual noise that occurs in nature: 

Vt{B} =B<^,t\^,t>B (4) 

Eq. (4) says that the statevectors with the largest norm are the most likely to 
occur. The total probability, according to (5), is 

1 f ~°° B 1 

V27rAt J- 00 

To see how these equations work, from Eq. (5) we note that the most probable 
B{tys occur if B{t) = 2Xtai or B{t) = 2Xta2 plus or minus a few standard 

deviations (At) 2. For small f. these regions overlap, but for large t they don't. 
For example, set B{t) = Bo{t) + 2Xtai, where the range of Bo{t) is a small integer 
x(Ai)2. Then, for large t, Eqs. (3), (5) become 
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\ip,t>B « .6|ai > e-3^^°W' + .8|a2 > e-^*[°i-«2l' 
Vt{B} « .366-5^^"^' + .64e-2^*["^-"=l' . 

Thus the probabihty associated with Bo{t) in this range approaches .36 and 
the statevector approaches ^ \ai >. Of course, a similar argument holds for 
B{t) = Bo{t) + 2\ta2, resulting in collapse to \a2 > with probability .64. For 
other ranges of B{t) the associated probability aproaches for large t. This, 
then, is the way B{t) acts as chooser: it chooses either 2\tai or 2At02 to fluctuate 
around, and this determines the collapse outcome. 

Collapse models such as this provide an explanation for the random results 
which occur in nature and arc unexplained by SQT: the result was this rather 
than that because the noise fluctuated this way rather than that way. What is 
left unexplained is at the next level: why did the noise fluctuate this way rather 
than that wayl A future theory may address this question by identifying the 
noise with something physical (e.g., gravitational fluctutions: see section 3.1) 
and having dynamics for it. 

It is useful to discuss various features of collapse dynamics in the context of 
this simple model before considering more sophisticated models. 



2.2a Density Matrix The density matrix describes the ensemble of statevec- 
tors which arise from all possible noises. For our simple example it is 



p{t) = -jl= r Vt{B}dB 
V In At J-oo 



\il;,t >B B< tp,t\ 

B< 1p,t\^,t >B 



= -^L=J dB\^,t>B B<^,t\ (6) 

= .36|oi >< ai| + .64|a2 >< 02] 

+.48 [|ai >< a2\ + |a2 >< e-^^"^""^)' (7) 

which follows from (4), (3). One clearly sees in the exponential decay of the off- 
diagonal density matrix elements that the collapse rate increases with increasing 
difference of ^'s eigenvalues. 

It should be emphasized that, while collapse dynamics gives rise to such 
density matrix behavior, it does not follow that such density matrix behavior 
implies collapse dynamics[42]. Thus, the unitarily evolving statevector 



= .6|ai > e-'^°W"i + .8|a2 > 6-'^"^"= 



where Bo{t) has probability density 
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Vt{Bo} 



e 



2\t 



has no collapse behavior, just increasingly random phases, yet at every instant 
its density matrix is equal to (7): 



It is often said that if the same density matrix arises from such different 
ensembles of statevectors as 

Ensemble A: \ai > (probability .36) plus \a2 > (probability .64) 
Ensemble B: .6|ai > e'^' + .S\a2 > e*^=^ (random 01,02) 

then the ensembles describe the same physics because they prcidict (statistically) 
identical experimental outcomes. The crucial point is that, for this argument to 
hold, it is necessary for the theory to say that experiments have outcomes. 

Suppose the state (1) represents the initial result of a measurement inter- 
action where \ai >, \a2 > describe macroscopically different apparatus states, 
and suppose A, B represent possible final apparatus ensembles. Each member 
of ensemble A represents a specific measurement outcome and each member of 
ensemble B does not. This conclusion is the same if a second apparatus is intro- 
duced to measure the state of the first apparatus and the two evolutions take 
place as did the first (collapse for A, noncoUapse for B). Thus one cannot make 
the argument cited above that ensembles A and B describe the same physics be- 
cause the argument depends upon both ensembles having a collapse evolution, 
which B does not have. 

In a recent stimulating series of papers, Mermin implies that an important 
lesson of quantum theory is that there can be no physical difference if there is 
no experimental difference [43]. This is in keeping with the lesson learned by the 
previous generation of physicists from the Ether issue. Well, each succeeding 
generation may well unlearn the lesson of the previous generation, e.g., when 
it's worth going to war. The position taken here is the direct opposite, that an 
important lesson of quantum theory is that just because we cannot measure a 
difference does not mean that there is no difference. After all, why should we 
be able to measure all that the world contains? Thus, consider the ensembles 
A and B but now take both ensembles as subject to the collapse dynamics (2), 
(4). Although the density matrix for both ensembles remains .36|ai >< Oi| -|- 
.64|a2 >< 02! during the evolution, the dynamics recognizes that the states that 
make up A and B are diff'erent and it treats them differently, leaving the states 
in A alone and collapsing the states in B. 
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2.2b Fourier Form It is useful and interesting to write the statevector evolu- 
tion (2) as a Fourier Transform: 

V27r J-oo 

This says that the collapse evolution can be viewed as a gaussian- weighted sum of 
unitary evolutions. This is useful because unitary evolutions are well understood- 
for example, one may use Feynman diagram techniques with them (see section 
4.2a). It also may be viewed more profoundly, as a possible natural generalization 
of the usual unitary evolution of SQT. In SQT the statevector not only represents 
the world but it carries the extra burden of representing a probability, so its norm 
must be one and its evolution is limited to being unitary. When the statevector 
represents reality and is freed from representing a probability it is also freed to 
have a more general evolution than unitary. 

The density matrix can likewise be written as a Fourier transform using (6) 
and (8): 

1 f°° 

p{t) = ^ drye-'''e-^''^V0)e*^^^. (9) 

J-oo 

Again, as discussed in the previous section, although this is a collapse evolution, 
the Fourier form of the density matrix displays it as a Gaussian-weighted sum 
of non-collapse (unitary) evolutions. 

2.2c Norm The statevector evolving according to (2) gets diminished in norm 
(e.g., from 1 to ~ .6 or f» .8 in the example). Further similar evolutions decrease 
the norm further. However, as we have mentioned, the norm has no importance 
in a theory with collapse (e.g., amplitude ratios, operator eigenstatcs arc norm- 
independent). Because the evolution (2) is linear in \tp >, the right side may 
be multiplied by a function F[t, B{t)] and the probability density (4) scaled by 
without physical effect. For example, if F = exp[B^(f)/4At], the evolution 
equation becomes 

so the final statevectors grow in norm (e.g., ~ expXtaf). Incidentally, this may 
be written as 

\ip,t >B= e Jo ^ <i* > 

showing that (2) has a time-translation invariant form. 

If, because of prior experience with SQT, one is more comfortable with nor- 
malized eigenstatcs, one may feel free to normalize the statevector at any time, 
or even set F =b< tp,t\ip,t so that the statevector is normalized at every 
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instant of time. However, then manifest scale invariance is lost and the evolution 
equation is nonlinear in \ip >. 

For evaluating collapse dynamics with a computer [44] it is useful to have 
the statcvcctor always of norm 1 since compiitcrs don't handle well numbers 
which get too small or large. Also computers most easily generate random noise 
functions which average around 0, and also most easily calculate dynamical 
evolution in small steps. To do this, one may consider evolution between t and 
t + dt and write the noise dB{t) in this time interval as 

dB{t) = dBo{t) + 2\dtA{t) 

^ ^ B<1p,t\A\'lP,t >B 
B<ll>,t\^,t>B 

thereby causing the probability density to be exp—dBQ{t)/2Xdt, belonging to 

an increment of Brownian motion. The statevector may be multiplied by an F 
chosen to normalize it. The result [45, 33] is the nonlinear equation 

d\(f),t>Bo={dBo{A-A)-X{A-Afdt}\(l),t>B ■ (10) 

In a recent papcr[46] the authors wrote that a linear evolution such as (2) 
is "merely a mathematical relation" and went on to say that "to be truly use- 
ful" states should be normalized. Computational usefulness should not be con- 
fused with physical importance. The evolution (2) with initial state (1) models a 
physical quantity dB/dt whose random fluctuations (both positive and negative) 
about one of the eigenvalues of A drive the statevector toward the associated 
eigenvector. In the evolution (10) with initial state (1), the noise fluctuates about 
zero with positive excursions enhancing the eigenstate with the larger eigenvalue 
and negative excursions enhancing the eigenstate with the smaller eigenvalue. 
Indeed, although (10) is computationally equivalent to (2) and (4), one might 
argue that (10) is "merely a mathematical relation" since a symmetrical efl:cct 
of positive and negative fluctuations about a mean value (such as is embodied 
in (2)) is the usual type of behavior for a random physical quantity. As another 
example of the questionable physical nature of such descriptions as (10), it has 
been shown [47] that such an equation can never lead to a relativistic collapse 
dynamics, unlike the situation with the linear equation (section 4). 

2.3 Less Simple Models 

The model upon which more sophisticated collapse dynamics is based has the 
statevector evolution 

\^P,t>^= Te-^/o '^*['"(*)-^^^(*)lV,0 > (11) 

(T is the time-ordering operator). This evolution is in what may be called the 
"collapse interaction picture," where A{t) = cxpiHtAexp—iHt is a Heisenberg 
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operator and the statevector only evolves because of the collapse dynamics. The 
probability density for w is 

Vt{w}=^<ij,t\iP,t>^ . (12) 

This is more complicated than the previous model because the statevector (2) 
is a function of B{t) but the statevector (11) is afunctional oiw(t). The probabil- 
ity that this noise lies between ■w{t) and w{t)+dw is Ylf,^Q{2TrX/dt)~'^dw{t')Vt{w}: 
one imagines the time interval t divided into infinitesimal segments dt, with each 
w{t') as an independent variable. 

However, the previous model is a special case of this, when = so A is 
time independent. Then, writing B{t) = jQdtw{t), (11) becomes 

This is just the evolution (2) multipled by a time-dependent norm factor. Apart 
from this factor, the evolution only depends upon B{t) and not upon B at earlier 
times so one may integrate the probability density over B at earlier times (using 
w{t') = B{t') - B{t' - dt)), obtaining the probability (4) of B{t) alone. 

One may think of (11) as describing collapse like (2) over each brief inter- 
val dt. The equation attempts collapse to A(t)'s instantaneous eigenstates, with 
a rate proportional to the squared difference of its eigenvalues. Since A{t) is 
changing with time, its eigenvalues provide a "moving target" for w{t) to fluctu- 
ate around. Collapse to one state occurs if the collapse rate characterized by A 
is faster than the transition rate between these states characterized by H. Most 
generally, there is a competition between collapse and Hamiltonian dynamics. 
But, often states do not appreciably change during collapse, so the behavior is 
well approximated by setting H = 0. 

2.3a Density Matrix The density matrix which follows from (11), (12) is 

= Te'^ lo *[^W®i-i®^(*)l'p(o) (14) 

{A Bp = ApB and T implies time ordering for operators to the left of p(0) 
and time-reversed ordering for operators to the right). 

2.3b Fourier Form The Fourier transform of (11) is found by regarding each 
w{t) as an independent variable with conjugate variable r}{t), resulting in 

Drje-^Io ^'''WTe^/o <i*'7(*)[-(*)-2AA(t)]|^^ ^ ^ ^^^^ 

-OO 
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expressed as a superposition of unitary transformations. 

The density matrix is easily expressed in Fourier form using (13), (15): 

Jo 
Jq 

(Tfj is the time- reversed ordering operator). Of course, if the 77's are integrated 
over, the result is (14). But, the form (16) displays the density matrix as a 
superposition of noncollapse unitary evolutions, each describing the interaction 
of the operator A{t) with Gaussian- weighted noise ri{t). 

2.3c Many A's To model the collapse of a system where states differ in a 
number of collapse-relevant features, each feature characterized by an operator 
An, the evolution (11) may be generalized to 

|V,*>.= Te-^Sn/;''*t-''«-2A^''Wl^|V,,0> . (17) 

The ^„'s are usually taken to commute at equal times. It is to their joint eigen- 
statcs that collapse tends, with a rate proportional to the sum of the squares of 
their eigenvalue differences. 

As an example of how this works, suppose the j4„(f)'s are time independent 
and the initial state is IV', >= X^jCi|i > where \i >= \af^ > |a-^^ >... and 
Anlal""^ >= a-"Vi"'' >• Then the analog of (3) follows from (17): 

|V,t >B= E^^l^ > e-^Sj^^W-^^*"'"'!^ . 

i 

Collapse takes place to |i > with probability |cip when Bn{t) « 2Ata-"^ Then 

the magnitude of the jth cigcnbasis vector decreases ^ exp —Xt [al"-* — Oj"'*]^, 
showing that each different feature in \i >, \j > increases the collapse rate. 

An equivalent form of (17) may be obtained in terms of new operators A'^{t) 
and noise functions w'^ (t) related to the old ones by a nonsingular matrix trans- 
formatfon An{t) = J2j Knj{t)A'j{t), Wn{t) = Y.j Knj{t)w'j{t): 

1^, t >^,= Te-^ Jo <^*K(*)-2AA;(t)]=G.,K.(t)-2AA^ Q ^ (^g^ 

where the matrix G = K^K is symmetric and positive. 

It is possible that one may choose K to be singular. For example, if A-^s 
eigenvalues arc all equal then Ai and wi play no role in collapse. Then one may 
choose Kij — 0, so G\j = and Ax and w\ are excised from (18). (Of course, 
there is one less w' to integrate over.) 
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2.3d NonMetrkovicin Evolution Another generalization [48] of the evolution 
(11) is 

\lP t >^= Te~^ ^ /^*dtidt2[«;(ti)-2AA(ti)]=G(ti-t2)[«'(t2)-2AA(t2)]2|^ ^ ^ ^^^^ 

This is a nonMarkovian evolution because it does not satisfy the usual Markov 
property that two evolutions, one from Iq to ti followed by one from ti to t, are 
equivalent to a single evolution from to to t. In Eq. (19), 

1 1"^ 

G{ti-t2) = — due''^^*'-*^^G{u) (20) 

where G(t — t') is real, symmetric under exchange of t and t' and positive, i.e., 
JJ^ dtdt'f{t)G{t - t')f{t') ~ |/(a;)|2G(a;) > for arbitrary f(t), so G{lo) must 
be real, symmetric and positive. When G{uj) = 1 so G{ti — ^2) = ^{ti — i2), (19) 
reduces to the Markovian (11). A reason for using the more general form (19) is 
to adjust the spectrum of w, which may be used to lower the amount of energy 
a system gets from w during collapse. 

To see how (19) gives rise to collapse, consider again the initial statevector (1), 
let A be independent of t and, for definiteness, choose G{t—t') = (a/2) exp —a\t— 
t'\ {G{lj) = o?l{J^ + o?)). Then (19) yields 



{B'{t)-2\a^{t-a-\\-e-'^')\Y 

n=l 

where S'(t)= / dti«;(ti)[l - -(e""*! + e-"(*-*i'] . 

For t » this becomes (3), whose description of statevector collapse was 
discussed in section 2.2. 

In the Markovian case, the probability rule (12) applied at any time T > t 
gives the same probability for the noise w{t). In the nonMarkovian case this is 
not so. In the present example 

2 

n=l 

For fixed t, as T increases past t, the probability of w{t) only "settles down" 
(becomes essentially independent of T) when T — t » . 

What does this mean? We interpret the probability as a measure of rational 
belief based upon present information [49, 50]. The future is not known at the 
present: present probability is conditional upon present time. Although the most 
we can ever know about {w} is given by Voo, at time t we use Vt given by (12) 
because it represents all that can be known at time t. In the above example, if 
is small compared to the collapse time, the difference is not significant. 
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2.3e Fourier Form of the NonMarkovian Evolution The Fourier transfor- 
mation representations of the statevector (19) and its associated density matrix 
are 

J —oo 

p{T) = /°° Dr]e~^^-f -fo '^*'^*'''^*^^o,Ti*'*">ii*") 

J — OO 

^ot(*1'^2) is the inverse of G over the interval (0,T): 

/ dtiG(i,ti)Go^(ii,i') = <5(i-t') (23) 
and depends upon the differences of T, 0, ii, t2- For our example, 

G,},{t,t')^{l-^^5{t-t') -L{5(T-t)5'{T-t') + 5{t)5'{t')} 

-^{5{T-t)5{T-t') + 5mt')}. 

The inverse of G over an infinite interval, = G~^{t — t'), is much easier to 
find than its inverse over a finite interval. Indeed, from (20) and (23) with the 
integral limits (—00,00), it follows that G^^{t — t') is the Fourier transform of 
1/G(a;). 

Because G is so simple, it is useful to note that Eqs. (21), (22) may be 
written with G^lp replaced by G^^ and the double integral limits replaced by 
(—00, 00). This can be seen from the identity 



-k J J^^dtdt'f(t)G{t-t')f{t') _ 

With f{t) nonvanishing only in the interval < t < T this becomes 



e-5/ j"^ dtdt'f{t)G{t-t')f{t') ^ 

which must be equivalent (manifest when the ry's are integrated over for time 
values outside the interval (0, T) ) to 
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/ dtdt'f{t)G{t-t')f{t') ^ 

A generalized form of (21), also a Gaussian- weighted superposition of unitary 
transformations, is[51] 

In the additional ?7-independent term, A interacts with the noise w through an 
arbitrary function F: its asymmetrical time integrals are required to make the 
trace of the associated density matrix = 1. This density matrix is readily shown 
to be 



p(T) — Te~^ -fo 'i*'i*'lMt)«>'^-^^Mt)]G'(t,t')[A{t')^i-mA{t')] 

g-iA J^^ dt[A(t)i81-l(8)A(t)l dt' F(t,t')[A{t')iS}l+10A{t')] ^^q^ 

{G'{t,t') = G{t-t') + dti J^' dhF{t,ti)G-'^{ti-h)F{t' ,t2)). This isaspecial 
case of an even more general nonMarkovian density matrix first discussed by 
Strunz[52], whose statevector evolution is based upon complex noise and allows 
nonHermitian operators. 

2.3f All-Time Dependent Evolution A variant of the nonMarkovian evo- 
lution (19) makes the statevector evolving from time to time t nonetheless 
dependent upon w for all time: 

lip t> =Te~^^ /^<iti<it2[j«(ti)-2AAo,T(ti)]G(ti-t2)[j«(t2)-2AAo,T(t2)]^^Q^ ^24) 

where ^o,T(ii) = ^(^i) for < < t and ^o,T(ii) = elsewhere. 

How is this to be interpreted? If the evolution is Markovian {G{ti — = 
5{tx — t^)) then the statevector (24) differs from the statevector (11) only by an 
inessential numerical factor 



~4X 

e 



i.e., although the expression for the statevector depends upon past and future 
w's, these have no effect. The interpretive issue lies with the probability rule (12). 
The probability for w{t) given by Vt{w) changes quite abruptly as T increases 



Collapse Models 15 



from less than t to greater than t. However, the interpretation of Vt{w) is 
no different than already stated for a nonMarkovian evolution, namely it is a 
conditional probability, conditioned upon present information. At unaccessed 
times Vt{w) just gives a "neutral probability estimate" of w{t), i.e., that of 
white noise with zero mean. 

If the evolution is nonMarkovian, there is a new wrinkle. The statevector 
(24) and its associated probability at time T depend upon in the future of 
T and the past of in a non- neutral way, at least for a brief interval (~ in 
our example) . There's nothing contradictory about this, it's just the way things 
are according to the model: at the statevector level the present depends at least 
a little bit on the future. 

However, the evolution (24) gives the density matrix 

p{T) = Dr}e~^^ ^ J^^dtidt2r,{ti)G-\tt-t2)ri{t2) 

J —CO 

.Te"'^'^ /o^<i*';(*)[^(*)®i-i'»^(*)]p(Q) (25) 

which is equal to the density matrix (22) (see explanation in previous section) 
which results from the evolution (19). Therefore, at the ensemble level, there is 
no effect of the future upon the present. 



2.3g Time- Smeared Evolution Another nonMarkovian evolution form equiv- 
alent to (24) is suggested by the observation that, since G is positive, we may 
construct G^: 

Gh{t-t') = — / da;e*"(*-*') [G{u))\ ^ 

so /_!°^(itiG^(f-ti)G^(ti-t') =<:;(*-*')■ For our example, [G(a;)]5 =a/(a;2 + 
a^)^ and G'^{t - t') = {a/TT)Ka[a{t - t')]. 

If w'{t) = J dtGi{t - ti)w{ti), As(t) = J dtGi{t - ti)Ao,T{ti) which are 
time-smeared variables are introduced into (24) we obtain 



Although this looks Markovian, of course it is not. Because T time-orders A{t), 
not As{t), two sequential evolutions are not equivalent to one evolution since, 
for t > t', [T As{t)][T As(t')] ^ [TA,{1:,)As{t')]. We remark that (26) may be 
useful in a relativistic model where spacelike smearing of an operator should be 
accompanied by timelike smearing. 
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2.3h Quantizing the Noise Although it may appear that the following is not 
a collapse model, I shall argue that it is. 

Consider an operator W{t) to represent quantized noise[41, 53, 54]: 



W{t) = dw [e-'^'a{w) + e'^'a^iu)] (27) 

whore a{u}), a^uJ) are annihilation and creation operators {[a{u)), a^{u!')] = d{w — 
w'), etc.) of a mode of frequency ui which can take on positive or negative values. 
This makes W{t) unusual for a quantum field, which usually only has positive 
frequencies. As a result, [W(t), M^(f')] = which also is unusual for a quantum 
field. But it is a nice property for a quantum field which has classical aspirations 
since a complete eigenbasis for W{t) is labelled by eigenvalues at each time: 



W{t)\w >=w{t)\w >, <w;|«;'>= ]J S[w{t) - w' {t)] . 

t=—oo 

If |0 > is the state with no noise mode excited {a{uj)\0 >= 0), one readily finds 

<«;|0>=e-^-/'-~'^*"'W. (28) 

dW{t)/dt commutes with W{t): one must look elsewhere for Ws conjugate 
operator. It is 

n{t) = -rT cLj \-e-^*a{uj) + e'^^a^u)] (29) 
SttA 

where [n{t), n(t')] = and [Wit), n{t')] = i6{t - t'). 
Now, consider the unitary statevector evolution 

\^,T> = Te-^2A/;dtAt)ii(*)|o > |V,,0 > 

= j Dw\w > Te"'^^-/ir'^*^W*?^ < > 1^,0 > 

= I Dw\w > Te-^/-~'^*['"(*)-''^'''^(*)^V,0 > (30) 

where the last step follows from (28). Eq. (30) says T >= / Dw\w > [ip, T >^, 
where \ip,t >w is identical to (24) (Markovian version). The standard probabil- 
ity measure of the state \w > li^^T >^ arising from (30) is equivalent to the 
probability rule (12). 

In section 1 I emphasized that such a superposition as (30) cannot represent 
a possible state of the world. Then, what is the meaning of (30)? 
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In his famous "cat paradox" paper [55], Schrodinger wrote "For each measure- 
ment one is required to to ascribe to the ■^-function (= the prediction catalog) a 
characteristic, quite sudden change. . . The abrupt change by measurement. . . is 
the most interesting part of the entire theory." (30) does not represent a pos- 
sible state of the world but, in Schrodinger's terminology, it is the "prediction 
catalog" of such states. Each \ip,t >^ in the superposition represents a possible 
state of the world: one of these is realized, the rest are not. Because of the dy- 
namics presented here, an unambiguous reality assignment to one of the terms 
in the superposition (30) is allowed which is not possible with SQT. Instead 
of measurement-induced (whatever that means) abrupt collapses (occurring at 
ambiguous times), collapses occur every dt sec. Each >„, "branches" into 
a family of possible \'i(),t -\- dt >w''s with the "label state" \w > distinguishing 
among them. The states \ip,t >^ in (30) describe the choices and the attacdied 
noise state |w > is the chooser. One may think of an Omar Khayyamish "moving 
finger" going from t tot + dt and choosing the actual |w >'s new w{t + dt) along 
the way. It is this feature which is either lacking or which is present in poorly 
defined form in interpretations of SQT without collapse dynamics. 

Eq. (30) has a number of features which facilitate this interpretation. 

1. For those enamored of environmental decoherence, the noise is effectively 
a universal environment (which does not have to be defined on an ad hoc 
basis for each physical situation). This is so desirable it suggests looking for 
a universal environment as the physical source for the noise (gravity is an 
attractive possibility — see section 3. Id). 

2. The states >w \w > do not interfere with one another. This allows 
the interpretation that one of them corresponds to reality and the rest to 
unrealized choices (one would think that an unrealized choice should not 
influence reality). 

3. The attached state Iw > acts as a complete label of the past history of the 
state \ip,t >yj, making identification of its evolution unambiguous. 

4. Each state Itpjt >w is physically sensible in that macroscopic objects are 
(almost always) localized (see section 3.1a). 

5. The states \w > t >w are always orthogonal because the \w >'s are always 
orthogonal. The states \ip,t >w representing reality are not always orthog- 
onal, e.g., if the initial state 11^,0 >= a\ai > +P\a2 > evolves continuously 
under different w's to \ai > or \a2 >, these states obviously cannot be or- 
thogonal during the evolution. (Indeed, it is my opinion that what has made 
a sensible SQT interpretation impossible is that the states which represent 
reality are not always orthogonal.) 

It is possible, in a future more sophisticated theory, that some of these attrac- 
tive features may be altered, with attendant complications. But, it is interesting 
that a collapse model based upon classical noise and a completely quantum model 
can be essentially equivalent. Although I shall not pursue this any further here, 
a similar completely quantum construction can be made for the nonrelativistic 
and relativistic collapse models in the next sections. 
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3 Nonrelativistic Collapse Models 

The first Galilean-invariant collapse model was constructed by Ghirardi, Rimini 
and Weber (GRW)[56, 57]. In their "Spontaneous Localization" (SL) model they 
introduced two important concepts which may be regarded as connected to the 
two parameters which characterize their model. 

The first parameter, A"-*^ w 10^^ sec « 3 x 10* yrs, accompanies the concept 
that collapse is caused by a physical process which acts on all particles. A is 
the collapse rate for a single particle, chosen small enough so that an individual 
particle in a superposition is negligibly affected, but large enough so that there 
is a dramatic effect on a macroscopic collection of n particles in a superposed 
state, namely collapse in (An)~^ sec. 

The second parameter, the mesoscopic distance a « 10~^cm, accompanies 
the concept that the collapse process narrows wavefunctions so that particles 
are spatially localized, but not localized too much. A particle's widely spaced 
wavefunction is narrowed to dimensions « n by collapse. Since a narrowed wave- 
function has increased energy by the uncertainty principle, if a were too small 
then particles would get too much energy from the process (see section 3.1c). If 
a were too large, then macroscopic objects would not be well localized by the 
process. 

I shall not discuss the SL model any further as it is not in one of the forms 
considered here. It has a problem (i.e., the process destroys the symmetry of 
the wavefunction under particle exchange) and has been superseded by the CSL 
model discussed in the next section. (See however references [58, 59] for reviews 
of this seminal model.) 

3.1 CSL 

In order to construct a Galilean-invariant model with features of SL in the form 
of a type we are considering, take the evolution equation (17) with many A„'s 
and choose n to be a spatial index so that An{t) — > A{x,t): 



Essentially, A{x,t) is proportional to the number of particlcs/vol in a spherical 
volume of radius a centered upon x — say, for now, that the particles are nucle- 
ons. In (31), N{z) = ^^z)^{z) is the nucleon number density operator, (^(z), 
^'''(z) are the nucleon annihilation and creation operators at z) and N{z,t) = 
expiHtN{z) exp—iHt is the associated Heisenberg operator. 
The statevector evolution 




(31) 



|^,T>^= Te 



dtdx[tu(x,t)-2AA(x,t)]^ 



(32) 
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and the probability rule (12) constitute what I called the "Continuous Sponta- 
neous Localization" (CSL) model[60, 45] (in appreciation of GRW's SL model). 
In CSL the classical field m;(x, t) chooses to fluctuate around the nucleon number 
density (x a constant, smeared over a volume ^ a^) of one of the states in a 
superposition, thereby causing collapse to that state. 

3.1a Collapse Rate in CSL It remains to explore consequences of the model. 
In doing so it is useful to work with the density matrix 

p(T) = Te"* r *'ix[A(x,f)®i-i»^(x,t)]^^^Q^ _ ^33^ 

Suppose that IV^O >= J^^, fifc(0)|nA- > where the normalized basis vectors 
|nfe > correspond to different nucleon number density distributions: N{pc)\nk >= 
nfc(x)|n/c >. Setting H = 0, the off-diagonal elements of (33) decay as 

< nj\p{T)\nk >= Cj(0)4(0)e-^ /"^"["^-("^-""("^l' (34) 

where afe(x) is given by (31) with N{z,t) replaced by nfc(z). 

Consider n particles in a small clump of size << a. Suppose the initial 
state is a superposition of such clumps with widely separated centers x/j , where 
|xj — Xfcl >> a. We may make the approximation nfe(z) nS{z — x^). Then 
(31), (34) give 

< nMT)\nk >« c,(0)cfc(0)*e-^^"' 

showing that the rate of collapse to one clump is proportional to n^. 

Now consider an extended homogeneous object of uniform number density 
n(z) = p, where the states \nj >, \nk > describe the object in two different 
places. Let their overlapping volume be Vo and their nonoverlapping volume be 
Vi. Then, by (31), aj(x) « {An)^a^p inside the volume of each. Outside the 
volume of each, aj{x) « 0. The integral over the overlapping volume Vo gives no 
contribution in (34), and we obtain 

< nMT)\nk >« c,(0)4(0)e-(4-)'^i^(-'p)(^iP) . 

Thus the collapse rate for such an extended object is proportional to the number 
Vip of nonoverlapping particles (and to the number of particles in a volume 
~ a^). Even if the object is small, a cube with sides 10~^ cm in length, with 
modest nucleon density p = 10^^ nucleons/cc it has a collapse rate of w 10~^ 
sec. This is the reason for the assertion that macroscopic objects are localized 
except for very brief intervals. 

But, what is meant by localized? A collapsed statevector, corresponding to 
an object considered to be localized here, generally has a "tail." This is a "small" 
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piece of the wavefunction which is not here and which is generally decreasing 
exponentially with time at the collapse rate (as evidenced by the above off- 
diagonal density matrix behavior). For a discussion of criteria for the "smallness" 
of the tails SCO [8, 9]. 

In CSL the collapse rate increases as superposed states differ more by having 
more particles in different places. In this way the interaction of a system with its 
environment may increase the collapse ratc[61, 62]. However, the prime initiator 
of the collapse is the system itself when it is in a superposition of spatially 
different states. 

Some examples of superpositions have boon given where the authors have 
felt that the superposed states are macroscopic enough so that collapse ought to 
occur (e.g., a superposition here+there of many photons[63]; a complex molecule 
in a superposition of excited plus unexcitcd states [64]) and criticized CSL for 
not predicting collapse in these cases. But it has been argued that there is 
no conflict with experience in these cases because collapse indeed occurs when 
the full experimental situation (e.g., observation of the photons by a human 
detector [65]; detection of the state of the complex molecule [66]) are taken into 
account: with CSL, opinion can be replaced by calculation. 

3.1b Nonlocality and Locality in CSL Collapse of an individual statcvector 
is a nonlocal process: something you do here affects something there. For exam- 
ple, if a particle is in a superposition of wavepackets here and there and you turn 
on a position measuring device located near here, then the wavepacket there is 
affected. Either the particle is detected so the wavepacket there disappears or 
the particle is not detected and the wavepacket there is all that survives. In 
either case the particle was neither here nor there before the measurement and 
is either here or there after the measurement, so the reality status of the particle 
there has been changed by an action here. 

The model clearly displays accord with Bell's theorem: to produce agreement 
with the predictions of quantum theory there is nonlocal influence. The nonlo- 
cality arises in two ways. One is from the direct product structure of states which 
requires the apparatus here to be multiplied by the particle state there so that 
collapse engendered by the former affects the latter. The other way is through 
the probability rule which requires high probability ■u;(x, i)'s everywhere. Thus 
a high probability w{x.,t) at the site of the apparatus is (nonlocally) correlated 
with a high probability ^(x, t) at the site of the particle. 

However, this nonlocality on the statevector level cannot be used to send 
signals from here to there because one has no control over the field w(x, t). This 
is evident by considering the density matrix in Fourier form: 



where Un{T) = Texp— i2A dtdx.r]{'x.,t)A{x,t). As is well known, such a uni- 
tary transformation does not allow nonlocal influences, nor does a superposition 
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of such. 

3.1c Experimental Tests of CSL Since CSL lias different dynamics than SQT 
it provides testably different predictions. A straightforward test might involve 

interference of a sufficiently large object[67, 68, 69]. For example, if an n-particle 
bound state wavefunction passes through a two-slit screen, the two emerging 
packets will fluctuate in amplitude (play the gambler's ruin game) according to 
CSL, so they will not have the same amplitude at a distance from the screen. The 
resulting interference pattern will therefore be different (diminished in contrast) 
from that predicted by SQT where the packets have the same amplitude. Such 
an experiment is difficult to do. 

A more readily performed experiment is to measure the increased energy of 
particles due to the collapse process[58, 70]. The expectation value of the energy, 
H{T) = TracciJp(T), can be found using (33). The potential energy part of 
H commutes with A(x, t) and, when the exponential in (33) is expanded in a 
power series in A the kinetic energy part of H gives / dx[A(x, t), [A(x, t),iJ]] = 
const X J dxN{x.) so all terms of higher order than the first vanish. The result 
for n particles of mass nip is 

iJ(r) = H(«, + 2Ar„^. 

This says 10^^ nucleons gain an average of fti.3 eV/sec which corresponds to a 
temperature increase of « .001°K. Although the average energy increase is quite 
small, infrequently a particle suddenly gains a large amount of energy and that 
can be looked for. 

We now note that different particles may have different collapse rates. To 
accomodate this we replace the particle number operator N{z,t) in (31) with 
'^^gaNa{z,t), where Na is the particle number operator for particles of type 
a: a particle's collapse rate is Xg^, with for electrons, gp for protons, gn for 
neutrons. 

Using (33), the probability /sec, F, of excitation of a bound state \tp > to an 
excited state |^ > (irrespective of the center of mass behavior) is found to be[71] 

r = ^\<ct>\Y,gaU)i^i-QM>f 

3 

to lowest order in (size of bound state/a)^. Here the sum is over all particles 
in the bound state, Xj is the position coordinate of the jth particle and Q = 
nfijyij/^- ruj is the center of mass operator. This rate vanishes identically 

if 5a(i) ~ rrij. 

A number of experiments currently being done for other purposes (e.g., dark 
matter searches, neutrino detection, proton decay, etc.) look for the sudden ap- 
pearance of energy in a volume of matter. To date, for purposes of testing CSL, 
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the most sensitive experiment has been one which looks for X-rays appearing in 
a slab of Germanium[72] . If a Is electron in a Ge atom is ionized by the collapse 
process, the result is an X-ray pulse of magnitude 11.1 keV (the Is electron's ion- 
ization energy, emitted as photons from the Ge ion decaying to its ground state) 
plus the kinetic energy of the ejected electron (converted to photons by the elec- 
tron's collision with Ge atoms). The excitation rate of such an electron is found 
from the equation above to have its largest value w 5000[((7e/<7p) — (m,, /rup)]'^ 
counts/keV-kg-day at ~ 11.1 keV, assuming that a has the GRW value and that 
A has the GRW value for nucleons (so gp = gn = 1)[73]. 

The best experimental limit [74] is that the rate at « 11.1 keV is less than 
w .2 counts/keV-kg-day, resulting in Qe/gp < ISme/wp. Thus, according to 
CSL and experiment, electrons collapse much less rapidly than nucleons. This 
is suggestive (but only suggestive, given the assumptions about the parameter 
values) of mass-proportional coupling, ~ m^/mp. This means that N{z,t) 
in (31) is to be replaced by M{z,t)/mp, where M is the mass density operator. 
Mass-proportional coupling suggests that collapse is connected to gravity. 



3. Id Collapse and Gravity It is hoped that collapse dynamics will some day 
be seen to arise in a natural way from another area of physics. It has frequently 
been suggested that this area is gravity[75, 76, 77, 78, 79, 80, 81, 35, 36, 37]. 
Failing gravity at present to move toward collapse, collapse may move toward 
gravity as follows. 

Let us make a change of variables of type illustrated in (18). We change from 
A{x,t) and w(x, i) to M(x,t) and w'{x,t):'^ 



mp(7ra^)4 J 

w{x,t) = /dzw;'(z,i)e-i^ 

mp(7ra^)4 J 

Now, Xa/c « 10~^^ is a dimensionless number not far from Gni^/hc « 10^'^*, as 
was pointed out by Diosi[77]. By substituting for the new variables in Eq. (32) 
and replacing A by Gm^/ha we obtain as evolution equation 



-i f dt<izdz'[«,'(z,t)-M(z,t)]Se 4a2 [„,'(z',t)-M(z',t)] 

\ip,T>^,= Te •'^0 \'W,To> 

(35) 

^ Unlike (18), this transformation is singular because G has some vanishing eigenval- 
ues: its eigenvalue spectrum is exp— o^k^/2 which vanishes at |k| = oo. This can 
be handled technically in various ways and, as discussed following (18), causes no 
essential difficulty. 
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We may think of (35) as the nonrelativistic hmit of a General Relativistic 
collapse formulation, where M is replaced by the trace of the energy-momentum- 
stress tensor (-^c^), and w' is regarded as the curvature scalar (xc^/SttG). That 
is, the collapse is due to fluctuations in the curvature scalar about its "classical 
value," an eigenvalue of the operator representing the traced stress tensor. 

Karolyhazy [75] was the first to suggest that fluctuations in the metric tensor 
could be related to collapse. Squires and I[71] gave a crude model for such fluctu- 
ations. In it, Planck mass {nipi) point particles randomly appear and disappear 
from the vacuum in the neighborhood of e.g., a proton, remaining on average for 
a Planck time with an average density of one "Planckon" per proton Compton 
volume {h/nipcf. It was shown that two quantities, the spectrum of fluctuations 
of the nonrelativistic limit of the curvature scalar (xG~^), w = (47rG')~^V^(/>, 
{(j) is the gravitational potential) and the energy increase of the proton (due 
to the random gravitational force —rripWcj) exerted by the Planckons), match 
the comparable CSL expressions with certain numerical coefficients, so that two 
equations were obtained for the two CSL parameters: 



A= r^^, a=[^\ « 1.4 X IQ-^cm. 



2(37r)i \tt^ J AnipC 



Although this model should not be taken very seriously, it does suggest that 
features of CSL might arise from an understanding of curvature fluctuations. 

Gravitional quantities other than curvature may be considered to give rise to 
collapse. Diosi[77] was the first to propose a gravitational collapse model of the 
form discussed here, with fluctuations in the scalar potential about its "classical 
value" causing collapse. His evolution equation may be written as 



(36) 

-. Mz,t), M„(z,t)= / dz' e-^^M{z',t) 

|x — z| J [2na-')2 

It was pointed out by Ghirardi, Grassi and Rimini[78] that Ma{z, t), not M{z, t), 
must be used in the expression for 0(x, t). If M(z, t) is used, a is effectively equal 
to the particle's Compton wavelength which means that particles get too much 
energy due to collapse. 

With a suitable change of variable, w{x,t) = G J dzw'{z,t)/\x — z|, (36) 
can look much like (35). Using Vx ■ Va' = — in the integrand of (36) and 
integration by parts one obtains 



I , ^ dtdxdx'[tu'(x,t)-M<.(x,t)]T— S_[u,'(x',t)-Ma(x',t)], , ^ 

which has a similar interpretation to (35). 
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Penrose's suggestion [76] that collapse arises from the gravitational energy 
associated with the difference of gravitational fields of superposed states may be 
implemented by 



, , „ ^ dt<ixdx'[w(x,t)-V<^(x,t)]^ , , ^ . . 

|^,T>^„=re ""■'■^o '^''''|V,To> (37) 

where fluctuations in the gravitational field w = g about its "classical value" 
causes collapse. 

The density matrices associated with the evolutions (35), (36) and (37) are 
respectively 



where /iM(z, t) = [M(z, i) 1 - 1 M(z, t)], and 



p(T) = Te~^ -fo d-td.^d,^'M{^'i) ^'''^''a^ ' ' '^'^^^''*V(0) (39) 
where A(j){z, t) = [0(z, i) 1 — 1 (g) (/)(z, f)], and 



p{T) = Te-^ r rf*'^-[^v<^(-.*rp(o) (40) 

where Z\ V(?!)(z, t) = [V0(z, t) 1 - 1 V0(z, t)] . Eqs. (39) and (40) are the same, 
as may be seen by integrating by parts in (39). For purposes of comparison with 
(38) they also may be written as 



p{T) =Te ^"^^ I"-" ' ■'o p(0) . (41) 

Consider the collapse rate (x?i) for a superposition of two states of differ- 
ent spread out » a mass distributions. According to (38), this is essentially 
^ G(ziM)^/a, where the sum is over cells of volume into which space is di- 
vided, and AM is the difference of the state's masses in each cell, i.e., the sum 
of the gravitational mass-diffcrcncc self-energy of each cell. According to (41), 
the collapse rate (x?i) is essentially the gravitational energy associated with the 
mass-difference (since the integral over s ~ 1). The collapse rate for (38) is, so 
to speak, a local gravitational energy and the collapse rate for (41) is a global 
gravitational energy. 

Anandan[81] has suggested that a natural generalization of Penrose's sug- 
gestion to General Relativity is that collapse arises because of differences in the 
connection FP^,. Unfortunately, [Z\P]^ is not positive definite so this suggestion 
cannot be implemented with the formalism presented here. 
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There have been other proposals regarding gravitational collapse. Collapse 
based upon energy differences {not energy density differences) was first proposed 
by Bedford and Wang[82] . It has received an elegant formulation by Hughston[36] 
which, however, is equivalent to (11) with A{t) = H and A^^ = ti^ipic? . The 
problem with such a proposal is that, in collapse as in real estate, location 
is everything and energy localization does not lead to spatial localization. The 
proposal of Pcrcival[35] is based upon energy localization for small distances: it's 
extension to large distances has not been completed. The proposal of Fivel[37] 
is based upon collapse acting to decrease a measure of entanglement and has 
also not been completed. One cannot say at present whether these two models 
localize satisfactorily. 

4 Relativistic Collapse Models 

One reason for trying to construct a special rclativistically invariant collapse 
model is that, like Mount Everest, the symmetry is there. But there are other 
reasons as well. It seems that collapse and relativity are intimately related. In 
nonrelativistic quantum theory, one cannot communicate to there by initiat- 
ing collapse /lere, for any experimental setup whatsoever. This is remarkable: 
why should a nonrelativistic theory be prevented from long-distance communi- 
cation in this way when one can communicate long-distance another way, by 
just sending a particle from here to there with sufficient velocity? It looks like 
this noncommunicability is a holdover from a relativistic theory where long dis- 
tance (superluminal) communication could make it difficult to have a consistent 
theory. 

It might be possible to have collapse in a preferred reference frame, i.e., the 

comoving frame of the universe, and yet have the experimental results (but not 
the theoretical structure) agree with relativity. Aharonov and Albert [83] have 
shown that an instantaneous collapse postulate in a preferred frame will not 
predict experimental results which disagree with special relativity. This assumes 
that collapse produces well defined and sensible states, which so far has not 
been achieved by any instantaneous collapse postulate but has been achieved 
by collapse models. However, collapse models also bring along with them other 
effects, and these may not be relativistically invariant in a preferred- frame col- 
lapse model. For example, the spontaneous excitation rate of atoms in CSL can 
be regarded as a clock whose ticking is a measure of the collapse rate.^ Then a 
moving atom's rate will be time-dilated and one may discern the special frame 
as the one in which the rate is fastest. 

Of course there are many phenomena which arc preferred frame dependent, 
e.g., the 2.1°K blackbody radiation. There is nothing intrinsically wrong with 
having collapse produce preferred frame phenomena. Still, it would be rather 
peculiar (although the universe is peculiar) for collapse, which allows quantum 

^ This was mentioned to me some time ago by Renata Grassi. It suggests that the 
collapse rate of any moving object is time dilated — see section 4.1d 
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theory to be compatible with relativity, to carry along with it phenomena which 
are not relativistically invariant. So, it is not unreasonable to guess that collapse 
may be completely relativistically invariant and to construct relativistic collapse 
models. 



4.1 RCSLl 

The first relativistic collapse model [79, 84] has a Markovian evolution: 



IV', cr >w= Te 



(42) 



(ctq, cr, represent spacelike hypersurfaces, x = (x, t), dx = dxdt and (j) is a. 
Heisenberg scalar relativistic quantum field). Here collapse occurs because of 
differences of scalar field amplitude (not because of differences of particle number 
density as in (32)). The reason for the scalar field is to obtain the parameter 
a in a natural relativistic way. a is the Compton wavelength corresponding to 
0's mass n = h/acfv leV. The Hamiltonian under which (p evolves includes the 
interaction g(j) -.■iptp: of a fermion field with (j) so that (j) "dresses" a fermion, i.e., 
surrounds it with an average Yukawa field of extension a. 

Collapse described by Eq. (42) works as follows. A Fermion in a superpo- 
sition of here plus there is entangled with its scalar field: \iI>,(to >= |Fermion 
/iere> I scalar here> + | Fermion i/iere> [scalar there>. w "chooses" to fluctuate 
around one of these scalar fields causing collapse to its state and, in so doing, 
it collapses the attached Fermion state as well. The collapse rate of a massive 
Fermion in a superposition of locations xi, X2 is [79]: 
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In the comparable CSL expression, ^g^a/lQ-K 
exp — (x — x')^/2a. 

This is the good news: collapse works well. 



A and exp — jx — x'j/a 



4.1a Vacuum Excitation The bad news is that collapse works too well. For 
simplicity, remove the Fermions from the model. The ^'-vacuum state may be 
written as a superposition of eigenstates of 4>{x). The evolution (42) acting on 
the vacuum proceeds to collapse it toward one of these eigenstates. This amounts 
to excitation of the vacuum: ^'-particles axe created. Each momentum k mode 
of the vacuum may be regarded as a harmonic oscillator in its ground state. It 
gets excited so that its average energy increases as H = ^iiXt, a modest rate of 
«1 eV/300 million years. However there are an infinite number of modes so the 
total average energy increases as 
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modes 



i.e., at an infinite energy /sec- vol rate. When Fermions are restored to the model 
they also are produced at an infinite energy/sec-vol rate. 

In a relativistic theory, when any vacuum excitation occurs it must be infinite 
excitation. (A particle of a particular four- momentum produced from the vacuum 
in one frame has a different momentum in another frame so, since all frames 
are equivalent, all frames must have particles of all momenta come out of the 
vacuum.) Therefore, to make a sensible relativistic theory all vacuum excitation 
must be eliminated. 

4.1b Pictures of Relativistic Collapse Before looking at the cause of vac- 
uum excitation and exploring how to eliminate it, I want to briefly mention a 
few features of the RCSLl collapse process. The elucidation of these and other 
features is mainly duo to Ghirardi and c;oworkers[84, 47, 85, 86] (see his con- 
tribution in this volume). The features mentioned here turn out to embody a 
proposal of Aharonov and Albert [87, 88] regarding instantaneous collapse in a 
(noninstantaneous) collapse model. 

Essentially, in each reference frame the collapse process works as in nonrela- 
tivistic CSL (except that the collapse time for a moving object is time-dilated — 
sec section 4. 2d). This means that the picture of collapse dynamics is frame 
dependent. Suppose a particle is in a superposition of wavepackets here and 
there and that apparatuses here and there are set up to detect the particle at 
the same time in a particular reference frame. In this frame at that time both 
apparatuses simultaneously trigger collapse: say, the packet here grows inside 
the apparatus here while the packet there fades away inside the apparatus there. 
However, in any other reference frame, only one of the apparatuses triggers the 
collapse and the packet not in it grows or fades before it reaches the other ap- 
paratus. Thus, unlike CSL, in RCSL what triggers collapse and when or where 
the collapse^ takes place; are frame-dependent. The model docs not give objective 
(observer-independent) answers to what, when and where but, of course, these 
are not measureable properties. 

4.2 Removing Vacuum Excitation in Lowest Order: the Tachyon 

To see the reason for the infinite vacuum excitation, and how to cure it, we begin 
by generalizing the evolution (42), using a combination of the generalizations (18) 



* For a more detailed presentation of the material in this and the remaining sections, 
see reference [89]. In what follows we set ?i = c = 1 unless mentioned otherwise. 



and (19): 



.4 
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I , _ -i^ f dxdx'[w{x)-2'y4>{x)]G{x-x')[w{x')-2'y(j>{x')]. 

\%p,a>=Te ^•'"o \ip,aQ> 



(43) 



where G{x - x') = — i-^ / dke'^<''-''">G{k'^) . (44) 
(27r)'* J 

The nonMarkovian evolution (43) reduces to the Markovian evolution (42) when 
G(fc2) = 1 (G(x - x') = 6{x - x')). 

The density matrix associated with (43), in Fourier form, is 

p(a) =T J £,^g-27/~^d^d^'r;(x)G-i(x-x')r,(x') 

-i27 [" dxrj(x)4>(x) i27 f dxr]{x)4>{x) 

■e ■'"o pycroje •'"o . (45) 

We see from (45) that the density matrix is the same as for a non-collapse unitary 
evolution where the scalar fie ld <j){x) interacts with a classical noise r]{x) with 
correlation function r]{x)ri{x') = {4:'y)~^G{x — x') and spectrum {A'y)~^G{k'^). 

As is well known, noise will cause a transition between two states if it supplies 
a four-momentum equal to their four-momenta difference. To create a particle 
out of the vacuum, the noise has to supply the four-momentum of the particle. 
For the Markovian evolution, G{k'^) = 1 so the noise r] supplies every four- 
momentum. To prevent r/ from creating ^-particles out of the vacuum one must 
set G(^2) ^ o_ 



4.2a Feynman Diagrams It is possible to use Feynman diagrams to visualize 
calculations because the density matrix (45) is expressed as the sum of unitary 
evolutions. Suppose at first that there are no Fermions. The perturbation ex- 
pansion of (45) provides two kinds of lines. One is I, representing a 0-particle 
of four-momentum k^ . The other is — representing the propagator GikF'): when 
the integrals over r] are performed, only even powers of t) give a contribution, and 
its correlation function G{x — x') acts like a propagator connecting two (/)-lines. 

The unusual feature is that the Feynman diagrams represent a density matrix. 
For example, the diagram showing the lowest order vacuum excitation of a single 
^particle is II. One must think of a vertical line down the middle of the diagram, 
separating it into the parts representing the terms to the left and right of the 
initial density matrix p(Tb)- To the left (right), the unitary transformation is 
time-ordered (time- reverse-ordered), so one must think of the four-momentum 
coming in from the future to the past (i) to the right and going from the past 
to the future (t) to the left. The only thing that can cross from right to left are 
G-propagator lines. 

In the diagram U, four- momentum is conserved because there is zero four- 
momentum in the initial vacuum state (below the lines) and zero four-momentum 
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going out: k"^ goes in at the right, then it passes through the G-propagator and 
goes out at the left. The density matrix associated with this diagram is 

p{T) = jG{„')T J ^^a\k)\0 X 0|a(k) 

(a(k) is the annihilation operator for the 0-particle). The resulting energy in- 
crease is 

Trace j (ikw(k)a1'(k)a(k)p(T) =7^(a*^)^^(^ j 

Thus, as mentioned earlier, if = the vacuum excitation vanishes. How- 

ever, this is only if there arc no Fcrmions, and of course wc must have Fermions 
since the whole point of collapse models is to collapse their states. 

4.2b Timelike Four-Momenta Cause Vacuum Excitation When there are 

Fcrmions, then Fcrmion pairs can be excited from the vacuum since Fcrmions 
are coupled to the (/>-particles. Pair production is represented by the diagram U 
with V V attached to the ends of the two <^lines. Since all four-momentum passes 
through the G-propagator. the diagram amplitude is proportional to G[(pi+p2)^] 
where pi , p2 are the Fermion pair four- momenta. To make this vanish we must 
have (5(fc^) vanish for k'^ > (2m)^. 

But, that's not the end of it. Pairs of (j) particles are produced via the diagram 
12 with ^ ^ on the ends of the 0- lines (q is a closed Fermion loop). To avoid this 
we must have G{k'^) = for fc^ > (2/^)^ = {2eVf. That leaves very little of the 
timelikc spectrum on which G{k'^) doesn't vanish. Indeed, it must vanish there 
too. When the Fermion is charged, photon production from the vacuum (via the 
previous diagram where a single photon emerges from the closed Fermion loop) 
occurs unless G(fc^) = for k'^ > 0. 

4.2c Enter the Tachyon By removing timelike four-momenta from the noise 
spectrum G{k'^), vacuum production of particles with timelike four-momenta 

is avoided. But, with only spacclikc four-momenta remaining in the spectrum, 
the immediate question is whether the nonMarkovian evolution (45) describes 
collapse and gives CSL in the nonrelativistic limit. The answer is yes. Timelike 
four-m,omenta excite the vacuum, and spacelike four-momenta cause collapse. 
Consider the c ^ oo limit of G: 



G{x-x',t~t') =7^7 / da;(ike'"(*-*')-'''-("-"')G(^4-k') 
(27r)* J 

-^b{t-t')-\. / dke-''^-("-"')G(-k2) . (46) 
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In this limit, G is Markovian. If we choose G{k^) to be nonvanishing at only a 
single spectrum four-momentum —fj,^: 



GiP) = 6{kl-k' (47) 
then it follows from (46) that 

(48) 

The nonrelativistic limit (48) is a fine replacement for the CSL Gaussian[90]. (For 
the record, the exact expression for G{x) is —{8TT^ax)~^Ni{x/a) for spaceHke 
X = (x^ — xl)^ and G{x) = —{4'K^ax)~^Ki{x/a) for timelike x = [x^ — x^)2, 
whose spacelike oscillations with period a and timelike decay with time constant 
a/c gives rise to (48)). 

(47) is the spectrum of a tachyon of mass n = 1/a. With this choice, the scalar 
field of mass /i is redundant: it isn't needed any longer to get the parameter a 
into the model. One may remove (p from the evolution equation (43), replacing 
it with N{x) =:^p{x)t/j{x): (then 7 is dimensionless) or the trace of the stress 
tensor M{x) = m ■.^j{x)^{x): (then 7 has the dimensionsion of G, i.e., mass~^). 
The vacuum pair production diagram V_V vanishes because the spacelike four- 
momentum of the tachyon propagator (formerly called the G-propagator) cannot 
equal the timelike four-momentum of the fermion pair. 

If asked to propose a natural nonlocal relativistic structure, it is likely you 
would mention the tachyon. It is interesting that the need here for the tachyonic 
structure did not arise from any abstract desire for it but rather from the needs 
to eliminate vacuum excitation in lowest order and to obtain the nonrelativistic 
CSL limit. Tachyons and collapse are two distinct theoretical structures requiring 
careful treatment lest superluminal communication and the difficulties of causal 
loops creep in: it is interesting that they are conjoined here. 



4. 2d Collapse Rate and Energy Production Rate We now have a model 
with evolution equation 



\ip,T >„= Te 



}- dxdx'[w(x)-'2-iN(x)]G{x-x')[w{x')-2'iN(x')], , „ , , 

4t I \ ) 1 y n y n \ 1 1 k '^\iP^To > (49) 



(we have used N rather than M to facilitate comparison with CSL's Eq. (32)) 
which is finite in lowest order, and we can calculate things. The associated density 
matrix in Fourier form is 



p(T) = T [ Dr]e~'^^-I^--^'^'"^'''''^''^'^ \x-x'Mx') 



Collapse Models 



31 



-i27 r dxJ7(x)[iV(x)(8il-l®JV(a:)] 
•e ■'^0 p{To) 



(50) 



PiTo) [ dxdx'G{x - a;')T[iV(a:), [N{x'), p{n)\\ + .. 



(51) 



CSL's comparable expressions are of the same form, where one replaces G{x — x') 
by 5{t — t') exp — (a; — x'Y j^o?, and N by the nonrelativistic particle number 
density operator. 

The Feynman diagrams corresponding to (51), describing collapse of a single 
Fermion in lowest order, are c| | and | \^ which give the rate of depiction of the 
initial state (coming from (51)'s N{x)N{x')p{Jq) and p{Tq)N{x)N(x') respec- 
tively) and |— I which gives the transition rate to the new state (coming from 



In the following calculations I have replaced N{x^ = ^p{x)^p{x) with N{x) = 
2rmlP'(x) where ip \s a, Boson field to simplify the calculation by avoiding Dirac 
algebra. The relativistic energy increase of the particle is found to be 



which is independent of the particle's four-momentum {dH{t)/dt is an invariant). 

For /i/m << 1 and 7 = Xp,~^ this is the same (apart from numerical factors) as 
the nonrelativistic rate SX/Sma^ of section 3.1c. 

The relativistic collapse rate of a particle in an initial state a\L > +0\R >, 
where \L > and \R > represent two widely separated identical wavcpackets 
with momentum-space wavefunction ^'(p), is characterized by the decay of the 
off-diagonal density matrix element: 



In the nonrelativistic limit one can set E = m in the integrand and the result 
is the same as the CSL result —a/?* A (see the equation following (34) in section 
3.1a). More generally, if the particle has a fairly well defined momentum we may 
make the approx imation |if'(p)p fa ^(p — Po), with the result that the collapse 



4.3 Finite RCSL Model for Free Particles 

Beyond lowest order there are diagrams with internal Fermion lines and, like 
the Hydra, vacuum excitation rears its head again. Pairs are produced from the 
vacuum via the second order diagram V^. 

Moreover, a free Fermion can gain infinite energy. Prom the first order dia- 
gram f-| (utilized in the last section), where the incoming and outgoing lines are 
on-shell, easy kinematics shows that an incident Fermion at rest gains a small 



N{x)p{To)N{x')). 




d_ 

dt 



< L\p{T)\R >= -al3*^^p^l + (£^)2 j dp||a^(p)|2 



rate, ~ y'l — Wg, is time dilated. 
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energy /x^/2m (« 10~^eV for an electron) and momentum Hs/l + {iJi/2mY by 
emitting a (negative energy) tachyon. However, a Fermion can gain arbitrary 
energy via the second order diagram P . 

There are two reasons for this. One is that the sum of four-momenta of 
two tachyons is capable of adding up to any four momentum, e.g., the sum or 
difference of two Fermion four-momenta. The other is that a Feynman Fermion 
propagator ^ (p^ — m? + ie)^^ is capable of carrying away from a vertex any 
four momentum, e.g., the sum or difference of a Fermion and tachyon four- 
momentum. 

At present I know of throe approaches to try solving this problem, only one 
of which I have spent some time at. In one (largely unexplored) approach, the 
idea is to limit the four-momentum a particle propagator can transfer by intro- 
ducing a rclativistically invariant cutoff, effectively "smearing" {:il){x)'ijj{x):^ 
J^^ dxidx2F{x,xi,X2) :'iij{xi)tp{x2)-)- A second (unexplored) approach has been 
suggested by Diosi (private communication). He notes that the density matrix 
for electrodynamics may be traced over the electromagnetic field, resulting in 
a more general nonMarkovian density matrix than we have been using here for 
RCSL[91]. However, it is a form for which Strunz[52] has given a collapse dy- 
namics (and for which a simpler collapse dynamics than Strunz's is given at the 
end of section 2.3e). This statevector and density matrix evolution describes just 
the particle behavior in electrodynamics (as in Wheeler and Feynman's classical 
action-at-a-distance electrodynamics) and is rclativistically invariant. The col- 
lapse dynamics replaces the decoherence due to radiation by particles. There is 
no energy creation problem — indeed, the extra F term effectively absorbs en- 
ergy, accounting for the radiative loss of electrodynamics. This is not a good 
candidate for a fundamental collapse model in that e.g., a massive body in its 
ground state in a superposition of here plus there will stay that way. But, it 
suggests that it may be possible to construct a satisfactory collapse model with 
the extra freedom of a more general nonMarkovian density matrix. 

4.3a Removing Time Ordering The approach I shall discuss is unconven- 
tional, but it may be simply stated. The evolution equation is (49) in Fourier 
form, except with the time ordering removed: 

— 27 dxdx'ri{x)G~^{x—x')r}{x') 

dxri{x)\w(x) — 2^N{x)\t , „ , , 

■c ^0 >^ ^ > ' ^ IV', To > • (52) 

Likewise, the associated density matrix is (50) with T removed. 

What does this achieve? The Feynman diagrams with or without time- 
ordering are the same. However, the off-shell propagator < 0\T'ip{x)'ip{x')\0 > is 
replaced by the on-shell propagator < 0\^p{x)^p{x')\0 >. That is, even when par- 
ticle lines are internal, each vertex |— describes a process where a real Fermion 
emits/absorbs a real tachyon, thereby gaining a limited amount of energy (it 
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is readily shown that AE/E < ji/m for jj, « m). Moreover, there is no pair 
production, real or virtual, at any vertex: \Z_ vanishes since a spacelike four- 
momentum cannot equal a timelike four-momentum: 

T 

r dx'G{x-x')i;l{x') — >0 as T ^ oo (53) 

{tp+{x) is the positive frequency — the annihilation — part of il>{x)). The result is 
a completely finite S-matrix: for large T, because of (53), 



p{T/2) = / dx,...dx2nJ2G{xi, - Xi,)...G{xi,^_, - Xi,J 

n=0 ^ •'-T/2 ^ 

•[^_(a;i)^+(xi),[...[V'-(x2„)V+(x2„),p(-T/2)]...]] (54) 

(C means the sum is over all combinations of pairs of indices; for simplicity 
I have assiimcd that no antiparticles are initially present and so have omitted 
terms). We see from (54) that there is no need for renormalization. One 
may also set N{x) — ^p-{x)ip+{x) in the evolution equation (52). 

Removing the time-ordering operation still gives an evolution which describes 
collapse. Indeed, in the usual test of collapse one sets H = and then, since 
operators are time-independent, time-ordering is irrelevant. 

It is worth remarking that removal of time ordering in a usual field theory 
results in the trivial S-matrix S = 1 since energy-momentum conservation pre- 
vents a real process at each vertex such as a charged particle emitting a photon. 
It is also worth remarking that reinstating time-ordering in e.g., (54) would re- 
sult in a nonrelativistic expression since Tip+{x)'tl>-{x') is not Lorentz invariant 
(■)/'+ (a:), \l)-{x') do not commute at spacelike separation). A non-trivial finite 
relativistic S-matrix without time-ordering is possible because of the tachyon. 

This model, however, is only for free particles. If the Fermion interacts with 
other particles, the free-particle field V should be replaced by the Heisenberg 
field '^ij. The removal of time-ordering only applies to the V^jf 's involved in the 
collapse interaction, not to those involved in e.g., the electromagnetic interaction, 
or else the usual electrodynamic physics will be destroyed. But this allows the 
Feynman off-shell propagator to connect two vertices, one where a tachyon is 
emitted, the other where a photon is emitted, with the result that the photon 
can get infinite energy. 



4.3b Time-Ordering and No Time-Ordering Removing time-ordering in- 
troduces nonlocality. How bad is it? 

To see what no time-ordering does, we begin by comparing the two pertur- 
bation series: 



T{e-^r''*^(*)}= v(-i)v^oT r dtr,.-- rdt2 rdt, 

Jo Jo Jo 
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g-iHo(T-t„)y^Q^ . . . g-i-f^o(t3-«2)y^Q^g-iHo(t2-ti)y^Q)g-ii?oti 

Jo Jo 
With time-ordering, particles evolve freely for time ti, exchange momenta, evolve 
freely for time t2 — ti, exchange momenta, etc., until time T. Without time- 
ordering, particles evolve freely for time ti, exchange momenta, evolve freely for 
time T — ti and then evolve backwards in time for time T, and then do it all 
over again. 

The hallmark of no time-ordering is backwards in time evolution. Repeated 
forwards and backwards in time evolutions can take a particle outside its light 
cone. However, as shall now be explained, in the present model this occurs with 
very small probabililty. 

4.3c Particle Evolution Consider the wavcfunction of a particle which starts 
out centered at x = and with momentum centered at zero as well. Consider 
how this wavefunction spreads in lowest order. This turns out to be the same 
as the spread of a classical ensemble of particles. In this ensemble, each particle 
starts out at x = with velocity v = 0. At random times in the interval (0, T), 
particles suddenly emit a negative energy tachyon in a random direction, thereby 
gaining velocity « {^/m)c (assuming /i << to). Thereafter each particle contin- 
ues moving freely forward in time for the rest of the time interval and, following 
this, it moves freely backward in time for T sec. 

A particle which received its impulse at t = travels a distance {p/m)cT 
by time T but, in its subsequent time-reversed evolution, returns to x = 0. At 
the other extreme, a particle which receives its impulse at time T ends up, after 
its time-reversed motion, at a distance {i.j,/m.)cT from the origin in the direction 
opposite to its velocity. Thus the particles end up distributed in a sphere of 
radius {^/m)cT, with velocities of magnitude {ii/m)c directed radially inward. 
These distributions describe the wavcfunction. (If the wavcfunction is that of a 
particle initially moving with some particular nonzero velocity, the description 
is the appropriately Lorentz invariant one, e.g., the spread of the wavefunction 
is Lorentz contracted in the direction of motion.) 

If the particle emits/absorbs n tachyons, it undergoes n successive evolutions 
as described above. The sum of ladder diagrams |— | + P +... essentially de- 
scribes the particle as undergoing a random walk with spread w {^i/rrijcT (this 
spreads faster than the usual random walk ~ Vt). The super luminal part of the 
wavefunction is negligibly small — a tail. 

4.3d Approximate Locality Consider two separate clumps of particles, L 
and R, at t = 0, and the light cones emanating from them. Assume the light 
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cones do not overlap at time t = T. Suppose the statevector evolution is unitary, 
\ip, T >= U {T)\ip, >. One may demonstrate locality, that anything going on in 
R during the interval (0,T) doesn't affect what goes on in L, by showing that 
TTaceR{U{T)p{0)W{T)} = U'L{T)TraceR{p{0)}U'^{T), where Z7^'s operators 
only act upon the particles in L. 

In a standard quantum field theory, where the operators in U are time- 
ordered, one may write U{T)\ip,0 >= J7^(T)f7£(T)|'0, >, where the operators 
in U'j^{T) {U'fi{T)) act only within R's (L's) lightcone, and the above result 
follows. One can do this because the operators within R's lightcone commute 
with the operators within L's lightcone and because the operators in U which 
are not in these lightconcs and which may not commute with R's or L's operators 
do not contribute to the evolution of \tp, T > (they act on the no-particle state in 
their region of spacetime and their contribution gets negligibly small for T >> 
?i/2mc^, the virtual pair production time). 

In the model presented here, the density matrix in Fourier form is a Gaussian- 
weighted superposition of unitary transformations ?7^, where the operators in 
U are not time-ordered. Because of the attendant forward and backward in 
time evolution, the operators in U which act on e.g., R's particles, occupy more 
than R's lightcone. Neglecting the tails, R's particles at time T are within R's 
lightcone, but the operators that got them there occupy a cylinder whose cross- 
section at any time is the cross-section of the lightcone at time T . While all 
the operators within R's cylinder do not commute with all the operators within 
L's cylinder, nonetheless the operators within R's cylinder do not act on L- 
particles, and vice versa. Therefore it is still true in this case that U (T^)|i/;, >= 
U'^{T)U'R{T)\ij;, >, where the operators in U'r{T) {U'^{T)) act only within R's 
(L's) cylinder, and the above-mentioned demonstration of locality (or rather 
approximate locality, since we have neglected the tails) goes through. 

This shows that it is possible to make a finite relativistic collapse model. 
As more such models are constructed, common features may be separable from 
model-specific ones. This can lead to a better understanding of what is necessary 
in such models and, perhaps, a model with particularly satisfying features may 
stand out. 

4.4 Concluding Remarks 

There is a big difference between a conditional statement and an absolute state- 
ment: "i/ you win the lottery then you will get ten million dollars" can't compare 
with "you have won the lottery and you get ten million dollars." 

The statements of SQT are conditional. Faced with the statevector ci|ai > 
+C2|(i2 >, SQT says "i/ this is the description of a completed measurement then 
the physical state is |ai > or \a2 >." But actually, what the if is conditioned 
upon, what the words "a completed measurement" mean, lies outside the theory's 
ken. SQT is not a complete description of nature because it fails to predict a 
physical phenomenon, namely that an event does — or does not — occur. 
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Phenomenological models are introduced into physics to describe phenom- 
ena that present theory fails to adequately treat. Collapse models are phe- 
nomenological models. Their statements are absolute. Faced with the statevector 

Cijoi > +C2|a2 >, the collapse model says that represents the physical state. If 
you wait a bit, it may happen that the statevector is unchanged, and that's 
that. Or, it may occur that the statevector rapidly evolves to \ai > or \a2 >, 
and again that's that. By "that's that" is meant, in all cases, that the statevec- 
tor represents the physical state: the model tells you whether or not an event 
occurred. 

A phenomenological model becomes convincing when it produces experi- 
mentally verified predictions beyond those engineered into its construction and 
achieves respectability when it is wedded to the rest of physics. So far, neither of 
these has happened for collapse models, but they are fairly young. However, the 
predictions arc there (section 3.1c) and perhaps, in the hints of gravity (sections 
3.1c,d) or tachyons (sections 4.2 — in string theory models, tachyonic modes of 
strings appear and arc gotten rid of) there may be the seeds of a connection 
to another physics domain. 
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